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Abstract 

We demonstrate the role of the mass asymmetry in the energy of vanishing 
flow by studying asymmetric reactions throughout the periodic table and 
over entire colliding geometry. Our results, which are almost independent of 
the system size and as well as of the colliding geometries indicate a sizable 
effect of the asymmetry of the reaction on the energy of vanishing flow. 
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1. Introduction 

In the search of nuclear equation of state as well as of nuclear interac- 
tions and forces, collective flow has been found to be of immense importance 
[3, 0, 0, S]- Among collective flow, transverse in-plane flow enjoys special 
status. During the last two decades, much emphasis has been put on the 
study of collective flow [H, 0, 0, Lots of experiments have been performed 
and number of theoretical attempts have also been employed to explain and 
understand these observations. As reported by |5( for the first time and later 
on by many others, collective flow is negative at low incident energies whereas 
it is positive at a reasonable higher incident energies. At a particular incident 
energy, however, a transition occurs. This transition energy is also known 
as energy of vanishing flow or balance energy. This energy of vanishing flow 
(EVF) has been subjected to intensive theoretical calculations using variety 
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of equations of state as well as nucleon-nucleon cross-sections [ll, 0, 0, 0, [Ef . 
This also includes the mass dependence of EVF which have been reproduced 
successfully by various theoretical models [H, H, 0| • 

Interestingly, most of these studies take symmetric or nearly symmet- 
ric reactions into account. Recently, FOPI group studied the flow for the 
asymmetric reaction of m Ca + 197 Au [6]. They noted that the flow in the 
asymmetric collisions is a key observable for investigating the reaction dy- 
namics. In contrast to symmetric collisions, where center-of-mass is one of 
the nucleus, this quantity is not known a priori in asymmetric nuclei ex- 
perimentally. Later on, FOPI conducted experiment on 58 Ni and 208 Pb 
In other class of studies, the total mass of the system was kept fixed as 96 
units whereas charge was varied Q. Theoretically, recently Kaur and Ku- 
mar [8| conducted a complete study of the multifragmentation by varying 
the asymmetry of the colliding nuclei. Asymmetry parameter rj is defined 
as (At — A p / At + Ap); where At and Ap are the masses of the target and 
projectile, respectively. All these attempts point towards a need for the sys- 
tematic study of the disappearance of flow for asymmetric colliding nuclei. 
Further as noted, asymmetry of a reaction plays dramatic role in heavy-ion 
collisions [9] . This happens because excitation energy in symmetric colliding 
nuclei leads to larger compression while asymmetric reactions lack the com- 
pression since large part of the excitation energy is in the form of thermal 
energy. 

Note that some isolated studies with asymmetric nuclei are already done 
in the literature where the reactions of 20 Ne+ 12 C, 20 Ne+ 27 Al, 20 Ne+ m Cu, 
58 Ni+ 12 C, M Zn+ 27 Al, 1 H+ 197 Au, 12 C+ 197 Au, 197 Au+ 12 C, 40 Ar + 207 Pb etc 
are taken into account 

0,0,1 

111 ]. Interestingly, none of the studies focus on 



the energy of vanishing flow in asymmetric colliding nuclei. To address this, 
we here present a systematic study of the EVF as a function of asymmetry 
of the colliding nuclei. While total mass of the reactions remain fixed, the 
asymmetry is varied by transferring the neutrons/protons from one nucleus 
to other. This study is performed with the quantum molecular dynamics 
(QMD) model which is explained briefly in section 2. 



2. The quantum molecular dynamics model 

The quantum molecular dynamics model (QMD) simulates the reaction 
on an event-by-event basis [2j. This is based on a molecular dynamics picture 
where nucleons interact via two- and three-body interactions. The explicit 
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two- and three-body interactions preserve the fluctuations and correlations 
which are important for TV-body phenomenon such as multifragmentation jij. 

In the model, the (successfully) initialized nuclei are boosted towards each 
other with proper center-of-mass velocity using relativistic kinematics. Here 
each nucleon a is represented by a Gaussian wave packet with a width of \f~L 
centered around the mean position r a (t) and mean momentum p Q (t) 0: 

Mr,P,t) = ^^ e H^(*)} a H c [«.(*W«I. (i) 
(27rL) 

The Wigner distribution of a system with A? + Ap nucleons is given by 

A T +Ap , 
f(r,p,t)= Y ^ e [-{r-Mt)}V2L} e [-{p~p a (t)} 2 2L/^] (2) 

with L = 1.08 fm 2 . 

The center of each Gaussian (in the coordinate and momentum space) is 
chosen by the Monte Carlo procedure. The momentum of nucleons (in each 
nucleus) is chosen between zero and local Fermi momentum [= ^2m a V a (r); 
V a (r) is the potential energy of nucleon a]. Naturally, one has to take care 
that the nuclei, thus generated, have right binding energy and proper root 
mean square radii. 

The centroid of each wave packet is propagated using the classical equa- 
tions of motion ji[ : 

df a dH 
dt dv n ' 



dH 



dt df, 



(3) 
(4) 



a 



where the Hamiltonian is given by 



,V2 



a 

Our total interaction potential V tot reads as j^] 

ytot ^/Loc _|_ yYuk _|_ yCoul _|_ yMDI 

with 

V Loc = trffa - r p ) + t 2 5{f a - f p )6{f a - f 7 ), (7) 
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V Yuk = t 3 e-^- r ^ /m / (\f a - T P \/m) , (8) 

with m = 1.5 fm and ts = -6.66 MeV. 

The static (local) Skyrme interaction can further be parametrized as: 

tLoc fp \ , a [P N 



Here a, (3 and 7 are the parameters that define equation of state. The mo- 
mentum dependent interaction is obtained by parameterizing the momentum 
dependence of the real part of the optical potential. The final form of the 
potential reads as 

jjMDI ^ tJn 2 [Wa _ + _ fp y (1Q) 

Here = 1.57 MeV and tg = 5 xlO~ 4 MeV~ 2 . A parameterized form of 
the local plus momentum dependent interaction (MDI) potential (at zero 
temperature) is given by 



The parameters a, (3, and 7 in above Eq. (11) must be readjusted in the 
presence of momentum dependent interactions so as to reproduce the ground 
state properties of the nuclear matter. The set of parameters corresponding 
to different equations of state can be found in Ref. 0]. 

3. Results and discussion 

For the present study, we simulated various reactions for 1000-5000 events 
in the incident energy range between 90 and 350 MeV/nucleon. In particular, 
we simulated the reactions of g 7 +ff Na (77 = 0.1), \ 4 N +^ Mg (77 = 
0.3), l°B +fi Si (77 = 0.5), and %Li S (77 = 0.7) for total mass = 40, 
f 8 Ar Ca (77 = 0.1), ™Si +| Cr (77 = 0.3), j° Ne +°° Ni (77 = 0.5), 

and l°B + 7 £ Ge (77 = 0.7) for total mass = 80, ™Ge Zr (77 = 0.1), 

54, 

26 1 



5 iFe +\™ Cd (77 = 0.3), f Ca +|° Te (77 = 0.5), and f 2 Mg +| 6 Ce (77 = 0.7) 



for total mass = 160, and IfCd +lf Ba (77 = 0.1), jgSV +Jf Dy (77 = 0.3), 
6 2 ° 8 Ni +}f W (77 = 0.5), and f 8 Ar P6 (77 = 0.7) for total mass = 240. 
The present study is for semi-central collisions (i.e. b/b max = 0.25). Note 
that in some cases, slight variation can be seen for charges. The charges are 
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chosen in a way so that colliding nuclei are stable nuclides. A soft equation of 
state with isotropic standard energy dependent cugnon cross-section (labeled 
as Soft* so ) and momentum dependent soft equation of state with standard 
energy dependent cugnon cross-section (labeled as SMD) are used in the 
present reactions. 

The energy of vanishing flow (EVF) is calculated using the directed trans- 
verse momentum < Pf r >, which is defined as: 

i 

where Y(i) and p x (i) are the rapidity distribution and transverse momentum 
of i th particle, respectively. 

In Fig. 1, we display the time evolution of the directed transverse flow 
< P* r >, for various asymmetric reactions with system mass (sysmass) = 
80 and 240 units. From the figure, it is clear that transverse flow saturates 
quite early during the interaction. A negative value during the initial stages 
signifies the dominance of the attractive interactions. One also sees well- 
known trends of transverse collective flow with increase in the incident energy. 
A clear transition from the negative flow to positive flow is also visible. From 
the figure, one can see that nearly symmetric reactions (rj = 0.1 and 0.3) 
respond strongly to the change with incident energy compared to highly 
asymmetric reactions. One notices that for r/=0.1 in the system mass 240, 
the change in the flow between 40 and 400 MeV/nucleon is 50 MeV, whereas 
for the same energy range, the flow varies between -2 MeV and +18 MeV 
(net 20 MeV) for r]=0.7. The cause behind is that with the increase in 
the asymmetry of a reaction lesser binary collisions take place resulting in 
lesser density and therefore, less response occurs for variation in the collective 
flow. On the other hand, for nearly symmetric reactions, the binary collisions 
increase linearly with incident energy, therefore, huge difference can be seen. 

It would be interesting to see how mass dependence character of the EVF 
behaves at a fixed asymmetry. In Fig. 2, we display the EVF verses combined 
system mass by keeping the asymmetry fixed. We notice that in all the cases, 
a perfect power law dependence (with power factor close to 1/3) can be seen 
for all asymmetries right from 0.1 to 0.7. All points lie on the line indicating 
that asymmetry plays a major role in the mass dependence. For a given 
mass e.g. system mass = 80 units, we see a variation of 50 MeV in EVF for 
rj varying between 0.1 and 0.7. This also explains why in the earlier mass 
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dependence studies though average behavior was a power law, individual 
EVFs were quite far from the average law [H,0]. There this happened because 
no control was made for the asymmetry of a colliding pair. Had all reactions 
been analyzed on a fixed asymmetry, one could have obtain perfect match 
with average power law. Further, a careful comparison of Soft and SMD 
reveals that the difference between EVF with both methods narrows down 
for light nuclei with larger asymmetry. In another calculations at b/b max = 
0.5 (not shown here), one noticed that MDI gives smaller EVF for lighter 
nuclei whereas EVF with MDI is enhanced for heavier nuclei. The difference 
is clearly visible for larger asymmetries. The cause of this behavior lies in 
the fact that MDI are more attractive compared to soft static equation of 
state at low densities whereas they turn more repulsive at higher incident 
energies. Based on this information, it is clear that for lighter nuclei at large 
asymmetry (and therefore, low densities) MDI suppresses the EVF. These 
observations are more visible at large asymmetries. 

In Fig. 3, we display EVF as a function of rj for a fixed mass equal to 40, 
80, 160, and 240 units. The results using the Soft 4SO (top panel) and SMD 
(bottom panel) EoS are displayed for clarity. In agreement with all previous 
studies, EVF decreases with increase in the mass of the system with both 
Soft ,so and SMD EoS. This decrease has been attributed to the increasing 
role of Coulomb forces in heavier colliding nuclei. As discussed in earlier 
figures, a sizable influence can be seen towards EVF with variation of the 
asymmetry of a reaction. For lighter masses, the effect of the variation of rj 
can result about 40 MeV change in the EVF. As noted in absolute terms, 
lighter nuclei are more affected compared to heavier ones. Overall, one sees 
that the effect of the asymmetry of a reaction is not at all negligible. It can 
have sizable effect which goes as power law with power factor close to 1/3. 

In Fig. 4, we display the percentage difference AEVF(%) defined as 
AEVF(%) = {(EVF^°-EVF rF=0 )/EVF rF=0 )xlOQ. Very interestingly, we 
see that the effect of the asymmetry variation is almost uniform throughout 
the periodic table. The mean variation is 0.7%, 4%, 13%, and 39% with 
Soft iso EoS and 0.3%, 2%, 6%, and 20% with MDI for rj = 0.1, 0.3, 0.5, and 
0.7, respectively. In other words, asymmetry of a reaction can play significant 
role in EVF and the deviation from the mean line can be eliminated if proper 
care is taken for the asymmetry of colliding nuclei. 

In Fig. 5, we display the EVF as a function of the impact parameter 
for different asymmetries. A well known trend i.e. increase in the EVF 
with impact parameter can be seen. Further, as demonstrated by many 
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authors, the impact parameter variation is less effected in the presence of 
MDI. The striking result is that the effect of mass asymmetry variation is 
almost independent of the impact parameter. 

4. Summary 

Using the quantum molecular dynamics model, we presented a detailed 
study of the balance energy with reference to mass asymmetry Almost 
independent of the system mass as well as impact parameter, a uniform 
effect of the mass asymmetry can be seen at the energy of vanishing flow. 
We find that for large asymmetries (rj = 0.7), the effect of asymmetry can 
be 15% with MDI and in the absence it can be 40%. This also explains the 
deviation in the individual EVF from the mean values as reported earlier. 
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Figure Captions 



FIG. 1. (Color Online) The time evolution of the directed transverse 
flow < P^ ir > as a function of the reaction time using different mass asym- 
metries. The left panel is for the reaction leading to sysmass = 80 units 
whereas right panel indicates sysmass = 240 units. The results for different 
mass asymmetries 77 = 0.1, 0.3, 0.5, and 0.7 are represented, respectively, by 
the solid, dashed, dotted, and dashed-dotted lines. 

FIG. 2. (Color Online) The EVF as a function of total mass (sysmass) of 
reacting partners. The upper panel is for Soft lso EoS whereas lower panel is 
for SMD. The results for different asymmetries rj = 0.1, 0.3, 0.5, and 0.7 are 
represented, respectively, by the solid squares, circles, triangles and inverted 
triangles. Lines are power law fit oc A T . 

FIG. 3. (Color Online) The EVF as a function of asymmetry parame- 
ter 77. The upper panel is for Soft* so EoS whereas lower panel is for SMD. 
The results for the systems having total mass A of 40, 80, 160, and 240 are 
represented, respectively, by the open squares, circles, triangles and inverted 
triangles for the case of SMD whereas shaded area in the case of Soft lso EoS 
represents the effect of charge variation leading to fixed value of charge and 
mass of the reaction. Lines are to guide the eye. 

FIG. 4. (Color Online) The percentage difference AEVF(%) as a func- 
tion of sysmass of reacting partners. The upper panel is for Soft* so EoS 
whereas lower panel is for SMD. The results of the percentage difference for 
different asymmetries 77 = 0.1, 0.3, 0.5, and 0.7 are represented, respectively, 
by the solid squares, circles, triangles and inverted triangles. Horizontal lines 
represent the mean value of AEVF(%) for each 77. 

FIG. 5. (Color Online) The EVF as a function of reduced impact pa- 
rameter for the sysmass 240. The upper panel is for Soft* 50 EoS whereas 
lower panel is for SMD. The results of the percentage difference for different 
asymmetries 77 = 0.1, 0.3, 0.5, and 0.7 are represented, respectively, by the 
solid squares, circles, triangles and inverted triangles. Lines are to guide the 
eye. 
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